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We construct a fully self-consistent Hartree-Fock-Bogoliubov theory that describes a spinless Fermi 
gas with long-range interaction. We apply this theory to a system of uniform dipolar fermionic 
polar molecules, which has attracted much attention recently, due to rapid experimental progress in 
achieving such systems. By calculating the anisotropic superfluid order parameter, and the critical 
temperature T c , we show that, "hign T c " superfluid can be achieved with a quite modest value of 
interaction strength for polar molecules. In addition, we also show that the presence of the Fock 
exchange interaction enhances superfluid pairing. 
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Introduction. — Recent experimental progress in ul- 
tracold polar molecules [lj has generated great interests 
in studying the properties and applications of such sys- 
tem. Applications associated with the internal energy 
levels of polar molecules range from quantum information 
processing Q to spin model engineering An equally 
intriguing direction is to focus on the external degrees 
of freedom the system of ultracold polar fermionic 
molecules with permanent electric dipoles represents an 
ideal setup to study dipolar effects in quantum degener- 
ate fermions [H, Q , as the dipolar interaction strength in 
these molecular systems is several orders of magnitude 
larger than that in atomic ones. 

Notably, two fundamental properties of the dipolar 
Fermi gases are superfluid pairing 0, Q and Fermi sur- 
face deformation [{J [n| , which are induced by the par- 
tially attractive nature of the dipolar interaction and 
the anisotropic Fock exchange interaction, respectively. 
Mathematically, the long-range interaction greatly com- 
plicates the calculation. As a result, pioneering works 
such as Ref. and Ref. 0, concentrated on each 
of these two features and also made further approxi- 
mations for simplicity. A quantitatively reliable fully 
self-consistent theory that includes both these features 
is lacking. 

The goal of the present work is to fill, at least on the 
mean-field level (which is believed to be reliable at low 
temperature for three-dimensional systems), this gap. In 
order to achieve this, we construct a self-consistent mean- 
field theory that takes full account of the interaction ef- 
fects. We show how this theory can be efficiently im- 
plemented by numerically calculating the superfluid or- 
der parameter and the critical temperature T c for su- 
perfluid transition. From our results, we show that ro- 
bust superfluid (with T c being a significant fraction of 
Fermi temperature) can be easily reached with ultracold 
polar molecules. We also investigated the interplay be- 
tween Fermi surface deformation and superfluid paring 



and show that the Fock exchange interaction enhances 
superfluid pairing via modifying the density of states. 

General theory. — We consider an ensemble of spin- 
less fermions with a general two-body interaction poten- 
tial U(r,r') = U(r',r) confined in an external trapping 
potential V(r). The second quantized Hamitonian reads 
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where ip is the fermion field operator, and \i is the chemi- 
cal potential. Denoting {f?(r)} as a complete set of single- 
particle eigenstates of —h 2 V 2 /(2m) + V(r) with eigenen- 
ergies £°, and the associated annihilation operator C v , 
Hamiltonian {!]) can be rewritten as 
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Performing the mean-field decoupling to the quartic op- 
erators, we obtain the effective mean-field Hamiltonian 

v v,v' 

-l[u h (v,v') + u f ( v ,v')}(clc v ,) - ^(WXC'tc't,) 

+±A( V M)Cfcl, + ±A*(r l ,ri')C n >C ri }, 
where the Hartree term Uh, the Fock term Uf and the 
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pairing term A are denned as 

Vl>V3 

H e g has a quadratic form and can therefore be diagonal- 
ized using the standard Bogoliubov transformation. 

Polar Fermi Molecules. — We now apply the general 
theory outlined above to a system of uniform dipolar 
Fermi molecules with dipole moment d = dz polarized 
along the z-axis. It is convenient to study this problem 
in momentum space. Instead of 77, we use the momentum 
k to label the single-particle states with e° k — ti 2 k 2 / (2m) . 
The interaction potential in momentum space is given by 



with 



U(q) = (4tt/3) d 2 (3 cos 2 q - 1) , 



(2) 



where q is the angle between q and the z-axis. 

From the symmetry of the system, at least for not too 
strong interaction strength, we anticipate that pairing 
only occurs between a particle with momentum k and 
another with momentum — k. In other words, the ground 
state has the usual BCS form: 

\gs) = Y\ C"k + v ^ C t C --k) |vacuum) . 

k 

Consistent with this ground state, the effective Hamilto- 
nian can be written as 



H cS = - [ Cl C_ k ] 
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where E Q = ±£[e(k) ~ U f (k)(CtC k ) - A(k)<c£cl k >], 
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A(k) = ^(7(k-k')(C_ k ,C k .) 
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Note that the Hartree term U h (k) = U"(0)X)(C£,Ck') 

k' 

vanishes as, for dipolar itneraction, U(Q) — 0. In ad- 
dition, it is easy to see that C//(k) = Uf(— k) and 
A(k) = -A(-k). 

The effective Hamiltonian (J3)) takes the diagonalizcd 
form 



Scff = [ £; ( k )7k7k - £(k) 7 _k7lk] + Eo 



in terms of the quasi-particle operators 
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where E(k) = e(k) 2 + |A(k)| 2 represents the quasi- 
particle dispersion relation. 

We remark that the quasi-particle dispersion E(k) may 
appear to have a similar form as that in the usual BCS 
theory for a two-component Fermi system with contact 
interaction. There is however a notable difference: In the 
usual BCS theory, the Hartree-Fock term is ignored as it 
represents a constant energy shift and can be absorbed 
into the definition of the chemical potential. By con- 
trast, here the Hartree-Fock contribution (for the uniform 
system considered here, only the Fock term survives) is 
anisotropic, due to the anisotropy of the dipolar interac- 
tion, and must be included explicitly. In fact, even for 
quite modest dipolar interaction strength, the Fock term 
has important effects and can lead to quite significant 
deformation of the Fermi surface 0, [l(| • 

At thermal equilibrium, we have (7^ 7k) = f(E(k)), 
(7-k7± k > = 1 - /(£(k)), where f(x) = 1/(1 + e?*) is 
the Fermi-Dirac distribution function. Consequently, the 
self energy term Q and the pairing term © take the 
following forms: 
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v ' ■ tanh - — 



2E(k') 



(8) 



It is known that the gap equation is ultraviolet divergent. 
The origin of the divergence can be attributed to the fact 
that the dipolar interaction potential used here [Eq. ^j] 
is not valid for large momentum. For large momentum, 
or equivalently for short distance, the dipolar interaction 
potential should be significantly modified due to repul- 
sion between electrons. Just as in the treatment of two- 
component Fermi gas with contact interaction, we need 
to regularize the interaction in the gap equation. This 
problem has been investigated by Baranov and cowork- 
ers In short, the bare dipolar interaction potential 
U(k — k') in the pairing term should be replaced by the 
vertex function 

r(k - k') = u(k -k')-J2 r(k - q)^jC/(q - k') , 



and the gap equaton should be renormalized as 

tanh 
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Equations ([8]) and ([9]), together with the number equa- 
tion N = ^2 k n(k) where 

n(k) = \u k \ 2 f(E(k)) + M 2 (l - f(E(k)) , 



FIG. 1: (color online) Contour of the momentum distribution 
function n(k) for Cdd = 1 at temperatures T = (a) and 0.3 
Tf (b). In (a), we draw the Fermi surface. The solid line is 
the Fermi surface obtained from the self-consistent calculation 
of this work, the dotted line is the one obtained from the 
variational approach developed in Refs. || [hJ. In (b), the 
lines from outside to inside correspond to n(k) = 0.1, 0.3, 
0.5, 0.7 and 0.9, respectively. 



is the momentum distribution function, comprise a com- 
plete description of the dipolar Fermi gas and need to be 
solved self-consistently. 

Results. — Now we present some results. First, let 
us consider a normal dipolar gas by taking A(k) = 0. 
Note that A = is always a solution to the gap equation 
©■ Fig.[T]ilustrates the momentum distribution function 

n(k) as a function of k z and k p - 



'k 2 + ky, for two 

different temperatures. Here the momentum is in units 
of the Fermi wave number of the non-interacting system 
kp = {Qir 2 !!,) 1 ^. The dipolar interaction strength is fixed 
at Cdd = 1 where Cdd = rndPn 1 / 3, /h 2 is the dimensionless 
dipolar strength Cdd = 1 corresponds to the RbK 

molecule created at the JILA experiment at a modest 



density of about 4 x 10 



At T = 0, n(k) 



0(— e(k)), where 0(.) is the step function. We draw in 
Fig-HUa) the contour of the Fermi surface. In Refs. [ilfioj]. 
we developed a variational approach and assume that the 
Fermi surface of the dipolar gas has an ellipsoidal shape: 

n(k) = (1 - a 2 k\ - k 2 p /a) . 

where a is the variational parameter characterizing the 
deformation of the Fermi surface. At Cdd = 1, we ob- 
tain a = 0.7769. In Fig. Q^a), the dotted line repre- 
sents the contour of the Fermi surface from this varia- 
tional calculation. As one can see, the variational re- 
sult matches with the full self-consistent calculation very 
well. At larger Cdd, small difference can be see between 
the two results. In general, the variational results exhibit 
slightly stronger deformation. The same conclusion has 
been reached by Ronen and Bohn 111 ], Fig. [ljb) shows 
the momentum distribution at T = 0.3 Tp, respectively, 
where Tp = Ep/ks is the Fermi temperature of the non- 
interacting system. At finite temperature, Fermi surface 
gets smeared out. However, the anisotropy of the mo- 
mentum distribution is still quite clear. 

Let us now turn to the discussion of the superfluid 
state. For simplicity, we take the first-order Born ap- 
proximation by replacing the vertex function T(k — k') 
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FIG. 2: (color online) (a) Gap parameter A(k) (in units of 
Ef) at T — for Cdd ~ 1. (b) The corresponding contour 
plot of the momentum distribution function n(k). The lines 
from outside to inside correspond to n(k) = 0.1, 0.3, 0.5, 0.7 
and 0.9, respectively. 
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FIG. 3: Ai(k) for Cdd = 1 at zero temperature. 



in the gap equation © by the bare dipolar interaction 
£/(k — k'). This should be a good approximation as long 
as the dipolar interaction strength is not too strong 0]. 
In Fig. [Ha), we plot the zero-temperature gap parame- 
ter A(k) for Cdd = 1- A(k) is an odd function of k and 
vanishes for k z = 0. As a consequence, the Fermi surface 
smears out except at k z — 0, as can be seen from the 
momentum distribution shown in Fig. d^b). The peak 
value of A reaches nearly 0.2Ep for this rather modest 
dipolar interaction strength, and occurs near k z = kp 
and k p = 0. To investigate the angular distribution of A, 
we note that: 

A(k) = A(fc,cos0 k ) = J2 Mk)Y w (cos9 k ), 

odd/ 

where k = |k| and due to the cylindrical symmetry of the 
system, only odd I and m — components are present. 
In Fig. El we plot A/(fe) for Cdd = 1- For small values of 
k (k < kp), A is dominated by the I = 1 (p-wave) com- 
ponent. For larger fc, contribution from higher partial 
waves may become important. 

Next, we illustrate the finite-temperature effects in 
Fig. [H The dashed line in Fig. |4ja) represents the chem- 
ical potential of the superfluid state as a function of tem- 
perature. It increases with temperature. In comparison, 
the chemical potential of the normal state (the solid line 
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FIG. 4: (a) Chemical potential (i as a function of tempera- 
ture for superfluid state (dashed line) and normal state (solid 
line). Energy and temperature are in units of Ef and Tf, 
respectively. Here Cdd = 1- The arrow indicates the location 
of T c . Superfluid state exists for < T < T c . (b) T c as a 
function of Cdd- Here the dots are numerical data and the 
smooth curve is a fit according to Eq. (|10[) . 



in Fig. \M, a )) is a monotonically decreasing function of 
temperature. Fig. 0Jb) shows how T c varies with Cdd- 
The solid line is a fit according to 



T C /T F = 0.8363 exp(-2.194/C dd ) 



(10) 



Some remarks are in order. First, from our calculation, 
we find that T c w O.lTp at Cdd = 1- If we were dealing 
with a two-component Fermi gas with contact interac- 
tion, such a critical temperature would correspond to a 
system inside the unitary regime. As we have mentioned, 
Cdd = 1 is a quite modest value for polar molecules. 
Therefore, typical polar molecules can easily reach the 
"strongly interacting" regime. Second, it is instructive 
to compare Eq. (TTT7|) to the critical temperature found by 
Baranov et al. [7] which in our notation takes the form: 



T C /T F w 1.44 exp 
= 1.44 exp( 
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4(6tt 2 ) 1 /3c (W 
.(-1.9887/C*,) . 



One can notice that the coefficients in the exponent agree 
quite well. Less agreement is found in the prefactor. This 
is, however, understandable as there are several differ- 
ences in our treatment. For example, Baranov et al. have 
included beyond-mean-field fluctuation and the contri- 
bution from the second-order Born approximation [l2j ]. 
while neglected the Fock term in their calculation. 

Finally, to reveal the interplay between pairing and 
Fermi surface deformation, we artifically turn off the Fock 
term in our calculation. We find that the presence of 
the Fock exchange interaction increases both the critical 
temperature and the magnitude of the order parameter 



by 20 ~ 25%. This enhancement can be understood in 
the following way. The presence of the Fock term causes 
an ellipsoidal deformation of the Fermi surface in such a 
way that it stretches the momentum distribution along 
the z-axis. As a result, the density of states near the 
Fermi surface is increased along z and reduced along the 
transverse directions. On the other hand, the dipolar- 
induced pairing is dominated by the p-wave symmetry, 
i.e., strongest in the z direction. Therefore, the Fock 
interaction-induced Fermi surface deformation tends to 
enhance superfluid pairing. 

Conclusion. - We have presented a fully self- 
consistent Hartree-Fock-Bogoliubov theory to study a 
system of spinless fermions with long-range interaction. 
We applied this theory to uniform polar Fermi molecules, 
calculated the superfluid order paramter and the critical 
temperature T c . Our work shows that: a typical Fermi 
gas of polar molecules can easily reach the "strongly in- 
teracting" regime with T c being a significant fraction of 
Fermi temperature Tp, and the Fock interaction has the 
effect of enhancing superfluid pairing. In the future, it 
will be of great interest to investigate the collective exci- 
tations of the superfluid dipolar Fermi gases, and the ef- 
fects of quantum fluctuations and the possibility of novel 
quantum phases that may arise at large dipolar interac- 
tion strength and/or in the presence of opitical lattice 
potential [1 31 ] - 
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